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We present a systematic description of the semiclassical kinetics of two-dimensional fermions in 
a smoothly varying inhomogeneous magnetic field B(r). We show that the nature of the transport 
depends crucially on both, the strength of the random component of B(r) and its mean value B. 
For vanishing B, the governing parameter is a = d/Ro, where d is the correlation length of disorder 
and Ro is the Larmor radius in the field Bo, the characteristic amplitude of the fluctuations of B(r). 
While for a <C 1 the conventional Drude theory applies, in the limit of strong disorder (a 3> 1) most 
particles drift adiabatically along closed contours and are localized within the adiabatic approxi- 
mation. The conductivity is then determined by the percolation of a special class of trajectories, 
the "snake states". The unbounded snake states percolate by scattering at the saddle points of 
B(r) where the adiabaticity of their motion breaks down. The external field B is also shown to 
suppress the stochastic diffusion by creating a percolation network of drifting cyclotron orbits. This 
kind of percolation is due only to the (exponentially weak) violation of the adiabaticity of the rapid 
cyclotron rotation in the field B, leading to an exponentially fast drop of the conductivity at large 
B. We argue that in the regime a ^> 1 the crossover between the snake-state percolation and the 
percolation of the drift orbits with increasing B is very sharp and has the character of a phase 
transition (localization of snake states) smeared exponentially weakly by non-adiabatic effects. The 
ac conductivity also reflects the dynamical properties of particles moving on the fractal percolation 
network. In particular, we demonstrate that the conductivity has a sharp kink at zero frequency 
and falls off exponentially at higher frequencies. We also discuss the nature of the quantum magne- 
tooscillations. We report detailed numerical studies of the transport in the field B(r): the results 
of the numerical simulations confirm the analytical findings. The shape of the magnetoresistivity at 
a ~ 1 found numerically is in good agreement with the experimental data in the fractional quantum 
Hall regime for the vicinity of half-filling of the lowest Landau level. 

PACS numbers: 71.10.Pm, 73.40.Hm, 73.50.Bk, 05.60,+w 



I. INTRODUCTION 

The transport properties of two-dimensional (2d) par- 
ticles moving in a spatially random magnetic field (RMF) 
B(r) oriented perpendicularly to the plane have attracted 
considerable interest in the last few years. This interest 
is largely motivated by the relevance of the problem to 
the composite-fermion (CF) description (l]J^] of a half- 
filled Landau level. Within this approach, the electron 
liquid in a strong magnetic field is mapped - by means of 
a Chern-Simons gauge transformation - to a fermion gas 
subject to a weak effective magnetic field. Precisely at 
half-filling, the average value of the Chern-Simons gauge 
field compensates the effect of the external magnetic field. 
The RMF appears in this model after taking static dis- 
order into account: fluctuations of the local filling factor 
induced by the random potential of impurities lead to a 
local mismatch between the gauge and external magnetic 
fields. A number of observations || of Fermi-surface fea- 
tures near half-filling give strong experimental support 
to the model of the effective magnetic field. Apart from 
the composite-particle models involving fictitious fields, 
2d electron systems with a real RMF can be directly 
realized in semiconductor heterostructures by attaching 
to the latter superconducting UH or ferromagnetic [HJ^] 
overlayers. 

The peculiarity of transport properties of 2d electrons 



in a random field B(r) shows up most distinctly in sys- 
tems with smooth inhomogeneities. The case of long- 
range disorder is most important also experimentally - 
since the compressible state in a half-filled Landau level is 
observed in high-mobility samples. In the latter, a large 
correlation radius of potential fluctuations, d, is deter- 
mined by a wide "spacer" between the electron gas and 
the doped layer containing ionized impurities. Likewise, 
inhomogeneities of the RMF created by the ferromag- 
netic overlayers in appear to be fairly long-ranged. 
The large value of the correlation radius d (as compared 
to the interelectron distance) allows to describe the elec- 
tron kinetics quasiclassically. 

It is well known that quantum interference effects may 
cause localization of noninteracting particles in an infi- 
nite 2d system even for arbitrarily weak disorder. This 
has been shown to be the case for charged particles in 
a random magnetic field [|||]. Specifically, the random 
magnetic field problem belongs to the unitary universal- 
ity class, with the localization length £ growing extremely 
fast with the dimensionless conductance g = cr xx / (e 2 / h) , 



£ oc exp(TrV). 



(1) 



These theoretical results are in full agreement with the 
recent extensive numerical study JlOj ] . According to (Q), 
already for g > 1.5 the localization length is larger than 
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any reasonable system size, and the quasiclassical ap- 
proach is fully justified. 

Let us stress that we consider the situation, in which 
the smooth random magnetic field constitutes the only 
type of disorder present. This should be contrasted with 
the starting point of Refs. |ll|], where the main contri- 
bution to the resistivity was assumed to be given by a 
short-range random scalar potential (treated within the 
relaxation time approximation for the collision integral in 
the Boltzmann equation), while a weak long-range cor- 
related random magnetic field was considered as a small 
perturbation. 

The purpose of this paper is to examine the transport 
in a long-range RMF in detail, with particular emphasis 
on the conductivity in an external (homogeneous) mag- 
netic field B = (B(r)) and/or at finite frequency ui. As 
we will see, the physics of the problem depends crucially 
on the values of B and uj, leading to a rich variety of 
different transport regimes. We complement the analyti- 
cal analysis by numerical simulations. The importance of 
the latter is due to the fact that in the most interesting 
part of the parameter space the transport is dominated 
by the phenomenon of percolation, so that only estimates 
"by order of magnitude" are available at the analytical 
level. 

In Sec. [n]we study the dc conductivity in a long-range 
RMF SB (r) at zero B. The character of the B = 
transport is determined by the parameter a = d/Ro, 
where Rq = vp(mc/ 'eBo) is the Larmor radius in the 
field Bo which is a characteristic amplitude of the fluc- 
tuations, vp the Fermi velocity. At a <C 1, the classical 
dynamics is of conventional diffusive nature, and the con- 
ductivity can be calculated in the Born approximation, 
which gives a xx ~ (e 2 nd/ wivf) / ot 2 , where n is the par- 
ticle density. With increasing a, a crossover to the per- 
colating kinetics takes place, and at a > 1 the conduc- 
tivity is determined by a small fraction of classical tra- 
jectories-so-called "snake states" [ |l2|]l3| . The extended 
snake states percolate through the strongly disordered 
system by winding around the lines of zero B(r) and 
yield a xx ~ {e 2 nd/ mvp) / a 1 / 2 L, which falls off slowly 
with increasing strength of disorder [here L — (In a) 1 / 4 
is a weakly varying logarithmic factor]. The crossover 
from the 1/a 2 to I /a 1 / 2 behavior of a xx at a ~ 1 is con- 
firmed by our numerical simulations. Furthermore, the 
latter allow us to find the numerical value of a xx for the 
CF problem (for which a ~ l/y/2 lies in the crossover 
region). 



In Sec. [II we consider the case of strong B. The trans- 
port through the snake states, examined in the previous 
section, is a dominant mechanism of the conductivity at 
strong disorder and zero mean B. Increasing B also leads 
to a suppression of the conventional diffusive motion and 
a transition to a percolation regime, even if a <C 1. The 
physics of this phenomenon is, however, quite distinct 
from the snake-state percolation. In the limit of large B, 



the trajectories are of the form of drifting cyclotron or- 
bits, and the dynamics is governed by an adiabatic invari- 
ant (magnetic flux through one cyclotron orbit). In the 
adiabatic approximation, the particles thus drift along 
the closed magnetic field contours and hence are local- 
ized. It is only the weak nonadiabatic scattering between 
drift trajectories that yields a finite conductivity. This lo- 
calization effect is similar to the formation of a "stochas- 
tic web" [ 14 la] in a slowly varying scalar random poten- 



tial in the presence of an external homogeneous magnetic 
field. The conductivity due to the weak nonadiabaticity 
falls off exponentially fast with increasing B. To calcu- 
late the function a xx (B) explicitly, one should consider 
the averaging procedure with care. An essential point is 
that the nonadiabatic transitions that determine the con- 
ductivity occur in rare fluctuations distributed sparsely 
along the percolating trajectories. These rare spots are 
characterized by anomalously sharp changes of the RMF. 
Using the optimum-fluctuation method, we find a Gaus- 
sian behavior of the conductivity in the limit of large 
B: In <j xx = — A(a) [B/Bq\ , where the coefficient A(a) 
scales as a 4 / 3 . Note that the conductivity in the CF 
problem (A ~ 1) falls off sharply beyond a small devia- 
tion from half-filling. It is worthwhile to notice also that 
the exponential falloff of a xx with increasing B is limited 
by inelastic scattering and by scattering on short-range 
static inhomogeneities - these scattering processes allow 
an electron to change drift trajectories and thus provide 
a mechanism of percolation competing with the nonadi- 
abatic transitions |l6|,f7[. 

While increasing B drives the system into the adia- 
batic regime at any a, the manner in which the conduc- 
tivity crosses over into this regime is qualitatively dif- 
ferent in the cases of weak and strong disorder. The 
situation is especially interesting at a > 1, where the 
transport regimes controlled by the snake states (weak B) 
and by the nonadiabatic dynamics of the cyclotron orbits 
(strong B) are separated by a sharp transition accompa- 
nied by an abrupt change of o~ xx . The transition occurs 
at B ~ Bq/ a 2 / 3 , where B becomes larger than the char- 
acteristic amplitude of the magnetic field at the nodes of 
the conducting network. At the critical point, the perco- 
lation network formed by the extended snake states falls 
apart into disconnected clusters, while the nonadiabatic 
scattering on the high-_B side of the transition is strongly 
suppressed and yields only a slight smearing of the criti- 
cal singularity. 

For a ~ 1, which is the regime relevant to the CF 
problem, we perform a numerical simulation to calcu- 
late the magnetoresistance p xx {B). The obtained curves 
are in good agreement with experimental findings: they 
show a weak positive magnetoresistance at low B, cross- 
ing over to a falloff of p xx with increasing B (in accor- 
dance with the exponential suppression of o~ xx in the adi- 
abatic regime). Further, we analyze the quantum oscilla- 
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tions of the resistivity and show that, in contrast to the 
conventional Shubnikov-de Haas effect in a short-range 
random potential, they start to develop only when the 
dimensionless conductivity <r xx / '(e 2 /h) drops down to a 
value of order of unity. 



In Sec. IV we discuss the transport in the RMF at fi- 
nite frequency u>. We find strong deviations of the ac 
conductivity from the Drude behavior, especially in the 
percolation regime, i.e. when a and/or B is large. At 
small u, we find a nonanalytical (oc contribution to 
the ac conductivity, which is determined by returns of the 
particle to the same spatial regions after a time ~ 1/M 
(and is analogous to the long-time "tails" found in the 
Lorentz gas many years ago ]l7j]). In the large- £? regime 
and at higher frequencies, the ac conductivity takes the 
form a xx (uj) oc |w| 3 / 7 , since 10 itself starts to determine 
the width of the percolating "stochastic web" responsi- 
ble for the conductivity. By contrast, in the snake-state 
percolation regime (a > 1, B = 0) a xx {u>) shows only a 
weak dispersion in the corresponding frequency range. At 
still larger frequencies the ac conductivity starts to drop 
exponentially reflecting the "ballistic" motion of drifting 
orbits (or snake states) on short scales. These analytical 
estimates are in agreement with our numerical simula- 
tions. 

Sec.0 summarizes our findings. The analytical results 
of Sections 0, |l| were partly presented in the Letter Q . 



II. DC TRANSPORT IN ZERO MEAN 
MAGNETIC FIELD 



A. Weak disorder 

We start with the simple case of a <C 1. In this limit, 
the CF trajectories are only slightly bent on the scale 
of d, so that the Born approximation is valid. Accord- 
ingly, for the transport scattering time one gets [p|Jig|l 
t^ 1 = Vp 1 (eB /mc) 2 J °° drf(r) = 2a 2 vp/d, where the 
CF effective mass m = fikp/vp is introduced. The Drude 
conductivity at zero B, a xx = e 2 nT tr /ra, then reads 



e 2 k F d 
'hJcT 2 ' 



a < 1 



(2) 



B. Dynamics of the snake states 

Let us now turn to the strong- RMF regime, a 1, 
keeping B = 0. The seemingly innocent assumption 
about the chaotic character of the particle dynamics, 
which enabled us to represent the conductivity in the 
form e 2 nTtr /in, is not valid anymore. Most particles are 
now out of play since they are caught in cyclotron or- 
bits drifting along the closed lines of constant B(r) ("van 
Alfven drift"). In the adiabatic limit, their drift trajecto- 
ries are periodic and so do not contribute to the conduc- 
tivity. Still, however large B is, there are classical paths 
which are not localized and percolate through the system 
by meandering around the lines of zero B(r). The con- 
ductivity is determined by the particles that move along 
these extended "snake states" Q (Fig. [I]). 



We begin by formulating the model to be studied. We 
consider non-interacting particles in the RMF B + 5B(r) 
with mean B and the correlator (SB(0)SB(r)) — B 2 f(r), 
where /(0) = 1. We assume that the function f(r) is 
characterized by a single spatial scale, which is the cor- 
relation length of the RMF. In particular, in the CF 
model with the electron density n equal to the charged 
impurity density rij we have B = (nc/e)(kp /v2d) and 
f(r) = (1 + r 2 /4<i 2 )~ 3 / 2 , where k 2 F = Aim (note that the 
electron gas is fully spin polarized near v = 1/2; we dis- 
card the spin degree of freedom throughout the paper). 
In this section we confine ourselves to the case of zero 
B. The RMF with zero mean is characterized by two 
length scales: d and the cyclotron radius Rq in the field 
Bq. Defining the parameter a = d/Ro, we can distin- 
guish the weak- RMF regime a«l, where the mean free 
path Z 3> i?o 3> rf, and the regime of strong fluctuations 
a 1, where one should expect drastic deviations from 
the Drude picture. We will explore these two limiting 
cases analytically. However, since the value of a corre- 
sponding to the CF problem lies in the crossover region, 
we will turn to numerical simulations in order to get a xx 
of the CFs. 




FIG. 1. Types of trajectories in a strong random magnetic 
field: drifting orbits along non-zero B contours and snake 
states near 5 = lines. Geometry of a snake state is char- 
acterized by the angle 8 (0 < 8 < n) at which the trajectory 
crosses the zero-field contour. Note that the direction of mo- 
tion of the snake state with 8 < 8 C (left) is opposite to that 
for 8 > 8 C (right), where 8 C ~ 131°. The width R s of the 
bundle of snake state trajectories is also indicated. 

Note that there is only one single percolating path on 
the manifold of the B(r) = contours; yet, the conduc- 
tivity is nonzero since the snake-state trajectories form 
a bundle of finite width, R s ~ d/a 1 ' 2 (see Fig. 0). The 
conducting network is made up of those snake states that 
can cross over from one critical zero-B line to another. 
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This coupling of two adjacent percolating clusters occurs 
near the critical saddle points of B(r), which are nodes 
of the transport network. The crucial role of the saddle 
points is that they break down the adiabaticity of the 
snake-state dynamics, as we are going to explain below. 

Everywhere except in small regions near the saddle 
points, the motion along the rapidly oscillating snake- 
state trajectories around the zero-B contours conserves 
the adiabatic invariant (see also Ref. [O]) 



I± = m (p ydy 



(3) 



Here y and y are the distance and the velocity in the 
direction perpendicular to the zero-B line, and the in- 
tegral is taken over one period of the oscillations. The 
conservation of the quantity I± can be established di- 
rectly by considering the evolution of the angle 9{x) the 
snake-state trajectory forms with the line of zero field 
(y = 0) at position x along this line as a consequence of 
the smoothly varying gradient b(x) = \dB(x,y)/dy\ y= o- 
The adiabatic invariant is parametrized as 



I ± (b,6) = Amvl /2 (2mc/eb)^ 2 F(6) , 



(4) 



where F{9) is a dimensionless function of order unity 
which can be found explicitly: 

F{9) = (1 - cos<9) f d£ v /(l-£ 4 )+cos(9(l- 1 e 2 ) 2 • (5) 
Jo 

Note that I±(b,9) may be written also as (e/c)$(fe, 0), 
where $ is the magnetic flux through the area encircled 
by the snake-state trajectory and the zero-B line in one 
oscillation period. We represent the equation dl±/dx = 
in the form of a scaling relation for the snake-state angle 



d9 
dhd) 



G{9) , G _1 (0) 



2±lnF(9) 



(6) 



This equation expresses the adiabatic invariance in terms 
of the fact that, given boundary conditions [9(x ) and 
b(xo) at some point xo], the angle 9 at a point x of the 
trajectory is completely determined by the gradient b(x) 
at the same point. Eq. (g) gives the asymptotic expres- 
sions for G(9): 



G{9) 
G(9) 



; 



9 -> 7T 



3(7r-0)ln^ 
Equations (|J), (0) tell us that 9{x) obeys the scaling 



fl(si) 
9{x 2 ) 



b{x x ) 
b(x 2 ) 



1/4 



in the limit of small harmonic oscillations 9 - 
singularity of G(9) in the opposite limit of 9 



(7) 
(8) 

(9) 



0. The 

■> 7r is a 



precursor of the bifurcation which accompanies the break 
away of the trajectory from the zero-f? line at 9 = w (see 
Fig. |, top). The functions F{9) and G{9) in the whole 
range of 9 are shown in Fig. |]a,b. 
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FIG. 2. "Serpentology" . Top: transformation of a snake 
state with large 6 into a drifting orbit with decreasing gradi- 
ent of the magnetic field; bottom: reflection of a snake state 
by a magnetic bottle-neck. 

The remarkable point to notice is that G{9) changes 
sign at some 9 = 9 C (which is ~ 131°). More specifi- 
cally, G(9) behaves singularly around 9 C , as (9 — C ) _1 , 
which corresponds to a maximum in F(9) at this point. 
This behavior of I±(9) means that the velocity of the 
snake states v s (9) (which is the average of x over one 
period) must change sign at 9 — 9 C , i.e., the snake state 
is "reflected" at the point x c defined by the equation 
9(x c ) — 9 C (Fig. ||, bottom). Indeed, as follows from Eq. 
(0), the constancy of I±[9(x)] cannot be maintained on 
both sides of the point x c . Note also that, in terms of the 
time evolution of 9, the change in sign of the function 
G(9) at 9 = 9 C means that the time derivative 9 retains 
the sign it had before the reflection. In fact, one can 
show, by solving the problem with constant gradient b 
exactly, that 



v s (9) - v F F'{9)- 



1 



\ cos 2 I 



sin0F(0) + cos 9F' (9) ' 



(10) 



i.e., v s (9) interpolates between v s (0) = vp and v s (tt) — 
—vp and vanishes at 9 = 9 C (see Fig. ||c). It is worth 
noting that the period of the oscillations T s (9) increases 
monotonically with growing 9: 



4 



Ts(0) 



„■> ! 2 -cx,t^) / 



(11) 



FIG. 3. The functions F(6) and G(0) determining the adi- 
abatic dynamics of the snake states according to Eqs. @-(^|), 
and the snake state velocity v s (9). 



i.e., T s (6) is equal to 2n(mc/ebvF) 1 ^ 2 at = and di- 
verges as 4(mc/e6i; J .) 1 / 2 ln[l/(7r - 0)] at -> tt. The 
"wavelength" of the snake states along the direction of 
propagation obviously reads Ax = \v s \T s , while the am- 
plitude of the oscillations in the perpendicular direction 
is given by Ay = 2vp(mc/ 'ebvp) 1 ^ 2 sin(0/2). 
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C. Snake-state percolation 

The adiabatic nature of the snake-state dynamics 
means that a typical trajectory is "trapped" between 
two return points x+ and x_ with 0(x±) — 9 C (Fig. |]a). 
Within the adiabatic picture, the drift motion in such 
a trap is periodic in time, as demonstrated in Fig. |]b. 
Hence, unless nonadiabatic corrections are taken into ac- 
count, these trajectories do not contribute to the dc con- 
ductivity. The nonadiabatic corrections for a typical tra- 
jectory with a slowly varying 0{x) are exponentially weak, 
so that the motion remains finite on an exponentially long 
time scale. Yet, this does not mean that a xx (a) is expo- 
nentially suppressed in the limit of large a. The point 
is that there are rare (but not exponentially rare) places 
along the zero-B contours where the adiabatic picture 
fails completely. These are regions where the contours 
pass near the saddle points of B(r), 
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FIG. 4. Snake state in a trap, (a) Real-space trajectory 
of a particle trapped between two bottle-necks. The scales 
of the x and y axes differ by a factor » 25: the figure is 
"compressed" in the x-direction. The dashed lines show the 
contours of the constant magnetic field, (b) Time evolution 
of the x coordinate. It is seen that the drift motion in the 
trap is almost periodic. 

B=0 




FIG. 5. Geometry of a saddle-point. 

Consider a snake state that is incident on a sad- 
dle point with the impact parameter p (Fig. 5). This 
means that the magnetic field at the saddle point is 
B sp ~ Bop/d and the distance Rmm at which the zero- 
B contour passes the saddle point is i? m in ~ y/dp. At 
the saddle point, there is an intersection of two lines 
of constant B(r) = B sp , while two zeio-B lines, along 
which the snake states can propagate, come within the 
distance 2i? m ; n from each other. Clearly, if i? m ; n is small 
enough, the snake state can change the zero-field con- 
tour. The angle 0, which characterizes the type of the 
snake-state trajectory, is then also changed, i.e., the adi- 
abaticity will be broken down upon "scattering" on the 
saddle point. To understand the parameters, consider 
first the case of p = ("direct hit"). The snake state 
propagates then along a straight line with decreasing gra- 
dient b(x) ~ Box/d, where x is measured from the saddle 
point. According to Eq. (H), 9{x) oc a; 1 / 4 decreases when 
the particle approaches the saddle point, while, as fol- 
lows from Eq. ([Tl]), the wavelength Ax diverges as x^ 1 ! 2 . 
The adiabatic picture is valid only as long as Ax(x) is 
much smaller than the scale on which the magnetic field 
changes, i.e., Ax(x) <C x near the saddle point, which 
gives the condition x 3> x c , where x c is the characteris- 
tic wavelength of the "last" oscillation before the particle 



hits the saddle point. For typical trajectories with ~ 1 
and Ax ~ (dRo) 1 ^ 2 at x ~ d, this condition fails at 
x ^ x c ~ d/a 1 / 3 . Now, as is clear from Fig. 6, whether 
the particle will be scattered to the left of the saddle 
point or to the right is determined by the initial phase of 
the snake-state oscillations. This sensitivity to the phase 
signals the breakdown of the adiabaticity. 
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FIG. 6. Scattering of a snake state at a saddle point. The 
particle may turn either left or right, depending on the initial 
conditions. 

We now turn to the case of finite p. At large enough 
p, the typical snake-state trajectory does not change 
the zeio-B line: the condition is that the angle 9(x c ), 
with which the trajectory comes to the saddle point, 
be much smaller than the ratio R m i n /x c . Substituting 
0{x c ) ~ (xdd) 1 ^ [see Eq. (§], we get p > d/a 5 / 6 . 
However, the adiabatic invariance is broken at the saddle 
point in a wider range of p: the condition for the curva- 
ture of the zero-B line to be large on the scale of the 
wavelength is R mm x c , which gives p < p s , where 

Ps - d/a 2 ' 3 . (12) 

Within this range, the angle 6 after the scattering, # ut) 
is typically of order unity even though the particle is in- 
cident on the saddle point with a small 9 (moreover, out 
depends strongly on the phase of the oscillations of the 
incoming trajectory). 

Now consider how the particles propagate between 
such non- adiabatic saddle points. The saddle points with 
p < p s are distributed sparsely along the zero-B trajec- 
tories with the linear density ~ p s /d 2 . Therefore, only 
a small fraction of the snake states can escape the adia- 
batic traps on their way between two such saddle points: 
most trajectories are localized in between. The snake 
state is not trapped if 9(x) < 9 C everywhere on its tra- 
jectory between the collisions with the saddle points. Ac- 
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cording to Eq. ([|), this is possible for trajectories with 
sufficiently small 9. Indeed, consider a snake state which 
has a small angle 9 <C 1 in a typical place with the gra- 
dient b ~ Bo/d. Typically, it will be able to travel a long 
distance, by far larger than d, until its angle reaches the 
value 9 C : this will occur in a fluctuation of the magnetic 
field with the anomalously high gradient b ~ Bo/d9 4 . 
Since the probability p(b) that the gradient at a given 
point exceeds some value b is determined by the Gaus- 
sian statistics, p(b) — exp(— b 2 / (b 2 )), we find that the 
state will typically run ballistically the distance L(9) 
obeying the equation L(8)p(B /d9 4 )/d ~ 1, which gives 
L{9) ~ dexp(0- 8 ). Hence, using L(d s )/d ~ d/p s ~ a 2 / 3 , 
we conclude that the states with the angles 9 <C 9 S , where 



-1/30-4 



> a- 1 / 12 



L(9 out ) ~rfx| exp^Jv* 2 ^)- 8 ] , 9°2 < a^ 2 



9 S ~ (In a) 



-1/8 



(13) 



will typically get through to the saddle point. 

We thus conclude that the particles with 9 <9 S prop- 
agate between saddle points "ballistically" with the lon- 
gitudinal velocity v s ~ vf, while others are simply out 
of play. Now we turn to construct the overall picture of 
the snake-state propagation. The scattering on a saddle 
point is actually a multi-step process. The fact that the 
angle 9 ou t is typically ~ 1 means that, having collided 
with a saddle point once, the particle almost inevitably 
returns back to it with a new angle of incidence 9': in 
effect, the trajectory "sticks" to the saddle point. How- 
ever, after many mappings 9 — > 9 out — + 8', the multi- 
ple reflections establish a stochastic distribution of the 
angle O ut characterizing the outgoing trajectory. Also, 
after many attempts the particle will go to the left or 
to the right with equal probability. This randomization 
of 8 ou t and of the direction of motion is clearly seen in 
the numerical simulation in p3]. Once the particle picks 
up the angle 6> out ~ 9 s /a 1 ^ , it will move ballistically 
until it reaches the next saddle point. Here, the factor 
0,-1/12 j g re i a ted to the fact that the angle 9{x) will in- 
crease ex a; 1 / 4 on the scale of d, so that the particle must 
have 9 ou t which is (xc/d) 1 / 4 times smaller than 9 S . At 
the new saddle point the whole process will repeat it- 
self. Since the saddle points are separated by the large 
distance ~ d 2 / p s , the average time it takes the parti- 
cle to move to the next saddle point is determined by 
the ballistic propagation between them, which requires 
the time % ~ d 2 j p s VF ~ (d/vp)a 2 ^ 3 , not by the mul- 
tiple attempts to "break away" with large # ut, which 
end in returns to the starting point. Indeed, assuming 
the full randomization of 9 out , we estimate the number 
of such attempts, until the particle picks up the angle 
$out ^ #out = a _1 ^ 12 (hio;) _1 / 8 necessary to reach the 
next saddle-point, as N ~ l/9^ t . According to what is 
said above, the initial condition 9 out allows the particle 
to advance the distance 



(14) 

Therefore, the average "waiting time" the particle spends 
in the unsuccessful attempts to reach the next saddle 
point is 

t w ~N [ d9 out L{9 out )/v F ~ (d/v F )a 2/3 /lna . (15) 



Thus, the total time it takes to get through from one 
saddle point to another is indeed determined by tb^> t w . 

D. Conductivity in a strong random magnetic field 

Now let us calculate the conductivity at a > 1. As 
was mentioned at the very beginning, most trajectories 
do not contribute to a xx since they follow periodic drift 
orbits. Next, we turned to consider a special class of the 
trajectories - the snake states. However, as we showed 
above, most snake states are also localized in the adi- 
abatic traps and only those with angles smaller than 9 S 
can propagate along the lines of zero B. At this point, we 
have to be concerned about the topology of the zeio-B 
contours. The first thing to notice is that all the contours 
are closed except one and this one percolating contour 
by itself cannot yield a finite conductivity. The conduc- 
tivity is nonetheless finite since the snake states in fact 
form a conducting network of finite width. The nodes of 
the network are critical saddle points, where two adja- 
cent percolating contours come close to each other. Note 
that most of the saddle points that the particle hits on 
its way between the critical ones only connect up small 
closed loops and so do not create a connected network. 
This happens only at the critical saddle points, where the 
snake states can cross over from one critical zero-B line of 
length L s ~ a 14 ' 9 d to another. We use here the results 
of the percolation theory (for a review see, e.g., fig] ): 
L s ~ d(d/ ' p s ) v+1 , where v — 4/3 is the critical exponent 
that controls the size of the critical cluster £ s ~ d(d/p s ) I/ , 
so that the ratio of the length of the trajectory and the 
distance from the starting point L s /£ s ~ d/p s . The char- 
acteristic distance between the nodes, i.e., the size of the 
elementary cell £ s , is then da 8 / 9 . On length scales longer 
than £ s , the particle dynamics can be viewed as fully 
stochastic. We estimate the macroscopic diffusion coeffi- 
cient as D ~ v s D s , where v s ~ L S R S 9 2 /^ 2 is the fraction 
of particles residing in the delocalized snake-states and 
D s ~ Q x vf/L s is their diffusion coefficient. Note that 
v s contains a factor 9 2 - since the density of the snake 
states is determined in the phase space parametrized by 
both the angle 9 and real-space coordinate: accordingly, 
one factor 9 S comes from the calculation of the fraction of 
the plane covered by the conducting snake states, while 
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the other describes their fraction in the 9 space. We thus 
have D ~ vpR s 9 2 and, correspondingly [BOfl, 



ki?d 



h oM*C 



r^in 1 / 4 , 



> 1 



(16) 



It is worth noting that the percolation enhances the 
conductivity: by comparison with the Born approxima- 
tion [Eq. (J2J) ] , the conductivity is ~ a 3 ^ 2 /C times larger 
(though the localization effects are strong and naively one 
might have expected the opposite) . Let us also note that 
Oxx given by Eq. ( |l6| ) is larger by a factor of ~ a 1 / 2 /C 
than that obtained for a ^ 1 in |21[] by using an "eikonal 
approach" . The fault in is not with the quasiclassical 
approximation itself, but with the method of disorder av- 
eraging, which neglects the localization of particles and 
the percolating character of the transport through the 
snake states. 

We now turn to the numerical simulation. To calcu- 
late the conductivity tensor components a^ u we evaluate 
numerically the classical current response function, 



' in/ 



e 2 pF 



di<« M (0)w„(t)) 



(17) 



where pp = m/2Tih is the density of states and the av- 
erage is taken over the disorder realizations and starting 
points of the trajectory. Typically, evaluation of the con- 
ductivities involved averaging over ~ 10 3 10 4 trajecto- 
ries. The numerical results for a xx in Fig. [7] fully confirm 
the analytical findings above. For small a, the results are 
in good agreement with the Born approximation formula, 
Eq. ML while at a ~ 1 a crossover to the a -1 / 2 behavior, 
Eq. (|l6|) , takes place. At a — l/y/2 (the value relevant 
to the CF problem at n = Uj and in the absence of im- 
purity correlations) we find a xx ~ 1.0(e 2 /h)kpd, which 
is a factor of ~ 2 larger than the Born approximation 
value. This improves the agreement with the experimen- 
tally found CF conductivity (defined as the inverse of 
the measured resistivity at v = 1/2), though the typi- 
cal experimental values of a xx are still larger than the 
one we obtain by a factor of ~ 2 — 3. This remaining 
discrepancy might be attributed to correlations in the 
distribution of the charged donors |22], which reduces 
the effective strength of the random potential and thus 
reduces a. The resistivity data in zero external magnetic 
field (as contrasted to zero effective magnetic field acting 
on CFs) indeed indicate that the model of statistically in- 
dependent impurity positions overestimates the amount 
of disorder |^,^2|. It is also worth noting here, in view 
of the controversy about the effective mass of the CFs 
], p^-|26| , that in the RMF model a xx at zero B [Eqs. 



(0),(16)] does not depend on m (neglecting the correc- 
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tions |p 7] related to the interaction between the CFs). 



FIG. 7. dc conductivity at B = 0, as a function of the pa- 
rameter a. The dashed and the full lines correspond to the 
theoretical asymptotics (^) and (|li|), respectively. Statistical 
errors do not exceed the symbol size. 



III. DC TRANSPORT IN NON-ZERO MEAN 
MAGNETIC FIELD 

We now consider the conductivity at finite B. Let us 
first discuss the case of a ~ 1, when the conductivity 
cr(B) can be parameterized as a function of the single 
variable B/Bq. As shown in the previous section, at 
small B/Bq we are at the crossover between the uncorre- 
cted diffusion and the snake-state percolation. Now, at 
B > 60 the particle dynamics is a slow van Alfven drift 
of the cyclotron orbits along the lines of constant SB (r). 
It follows that the conductivity is determined by a per- 
colation network of trajectories close to the SB(r) = 
lines. From the point of view of topology of the net- 
work, the situation is thus similar to that at zero B and 
a> 1. What is crucially different, however, is the mech- 
anism of the percolation. Specifically, at large B there 
is no stochastic mixing at the nodes of the percolation 
network: unlike the snake states at B — 0, the rapidly 
rotating cyclotron orbits pass unharmed through the crit- 
ical saddle points of SB (r) without crossing over to the 
adjacent cell. In the high-i? limit, the mixing occurs on 
the links of the network and is only due to the weak 
scattering between the drift trajectories. 

In order to calculate the conductivity at B ^> Bq, we 
should integrate out the fast cyclotron rotation, taking 
care not to lose the effect of the nonadiabatic mixing. 
Specifically, we have to go beyond the standard sepa- 
ration of the fast and slow degrees of freedom, known 
as the drift approximation. The parameter that governs 
this separation is 5/d, where S is a characteristic shift of 
the guiding center after one cyclotron revolution. The 
drift approximation is represented as a power series in 
5/d < 1. In our problem at a ~ 1, this parameter is 
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the ratio Bq/B. Therefore, if B 3> Bq, the adiabatic de- 
scription is good on microscopic scales. The key point, 
however, is that the conductivity is strictly zero at the 
level of the drift approximation - since the drift orbits 
are periodic in the thermodynamic limit. The effects 
that break the adiabatic invariance and lead to the tran- 
sitions between the drift orbits are exponentially weak at 
<J/d«l. 

A. Single-impurity scattering 

The problem of the scattering between the drift trajec- 
tories in the static RMF Jl8), as well as a similar prob- 
lem for a random scalar potential, considered recently in 
fl4|| , is a particular example of the broad class of prob- 
lems dealing with non-conservation of an adiabatic in- 
variant. Despite the general interest of this problem, 
any systematic expansion capable of giving the scatter- 
ing rate beyond the exponential accuracy, has proven to 
be a tough exercise. To consider a transparent example, 
we formulated and solved parametrically exactly a single- 
scattering problem |^9| . Specifically, we introduce a weak 
homogeneous gradient e of the background magnetic field 
and consider the interaction with an "impurity" modeled 
by a spatially localized perturbation 5B(r) of size d, so 
that the total field 



B(r) =B[l + e(y/R c )]+6B(v) , 



(18) 



where R c is the cyclotron radius in the field B. The 
guiding center coordinate y averaged over the cyclotron 
orbit, p = (y) c , plays the role of an impact parame- 
ter. The particle entering the system at x — — oo with 
(y) c = Pi wu ^ leave it at x = oo along the trajectory with 
(y) c = Pi + A/0, where Ap is the nonadiabatic shift we 
are interested in. In this single-impurity scattering prob- 
lem, the shift is a perfectly well-defined quantity. To first 
order in SB, the exact solution is given by 

t <55[r (i)l 



Ap = 7/ 



dt- 



1 = LJ± Vo(t) 



(19) 



Here ro(i) is the unperturbed trajectory for SB = 0, 

7 (e) = uo c (y (t) ^ dt'x Q (t')/y 2 (t')^ , (20) 

lj c = eB/mc, and the angular brackets denote averaging 
over one cyclotron period. In the limit e — -> the constant 
7 — > — 1. We can further simplify the model by assuming 
that SB is a function of x only. In this case, the integral 
in Eq. ( |l9| ) is evaluated at e <C 1 by the saddle-point 
method to give 



Ap 



2v F 



e I 2 

cos 

B V 7T \ e 



COS (fo 



(21) 



00 



I dtcosuj c t SB (—VFt 
-00 ^ 2 



-R, 



This formula expresses the nonadiabatic shift in terms of 
the asymptotics of the Fourier transform of the smooth 
function 5B{x) - thus demonstrating explicitly the ex- 
ponential smallncss of Ap. It shows that the parameter 
which governs the exponential falloff of Ap is d/S ^> 1, 
where 8 = ireR c , while the ratio d/R c may be arbitrary. 
Note that the pre-exponential factor happens to oscillate 
wildly as e — *■ 0. These oscillations are geometric reso- 
nances due to the commensurability of two length scales 
R c and S. Remarkably, the series of the resonances is 
defined by the properties of the unperturbed solution 
( "self-commensurability" ) and not by the shape of the 
scatterer. This means that the oscillations are damped 
with increasing strength of the perturbation |2^] - since 
the resonance condition cannot be met simultaneously 
everywhere on a strongly perturbed trajectory. Another 
peculiar feature of the nonadiabatic shift is its sensitivity 
to the phase ipo of the cyclotron rotation of the incident 
electron (Ap oc coscpo)- 

In the CF problem, a charged impurity located at 
a distance d from the plane occupied by the elec- 
tron gas creates the axially symmetric perturbation 
SB(r) — SB d 3 (r 2 + tf 2 )~ 3 / 2 . Because of the branch 
points at r = ±id in this expression, the integrand 
in Eq. ( |l9| ) will contain the exponentially small factor 

/ 2ui c 



suit reduces to 



Pi] 2 + d 2 I . The lengthy general re- 



5B VRcd 
Ap ~ snd -= — 

B 8 

/2 7T 

x cos ipo cos I — 1 exp 



(22) 



2ird 



at pi = in the limit d R c . Equation (|22j) reflects 
the features of the non-adiabatic shift discussed above: 
the exponential smallncss, the oscillations with chang- 
ing e, and the oscillatory dependence on the phase <po- 
These results were confirmed by numerical simulations 
in Ref. (29). Note that Eq. @ implies that the drift 
trajectory is only slightly perturbed by SB(r). 



B. Optimum fluctuation 



In the transport problem one has to average over an 
ensemble of impurities. What is crucial for the averag- 
ing process is that the drift velocity is itself determined 
by the fluctuations of the impurity field. The non-linear 
problem gets therefore much more involved as compared 
to the single-scattering model above, but the principal 
features of the nonadiabatic scattering remain unchanged 
and the main message can be simply stated: Because 
of the exponentially strong dependence of the shift on 



9 



the parameters of the single scatterer, the conductivity 
is determined by rare places with an anomalously high 
rate of nonadiabatic transitions jl4|,[l8| . Accordingly, one 
can neglect correlations between consecutive transitions 
from one drift orbit to another and each nonadiabatic 
shift can be considered independently. Since the nonadi- 
abatic scattering rate increases as the drift motion gets 
faster, the effective scatterers, sparsely distributed along 
the percolating trajectories, are characterized by anoma- 
lously sharp changes of the RMF. The problem now is to 
find the density and the parameters of these scatterers. 

The nonadiabatic shift Ap = Ap x + iAp y (in complex 
notation) after one scattering reads 



Ap = v F / dt^ t+itp0 A{t) 



(23) 



where the smooth function A(t) varies slowly on the scale 
of lu^ 1 and is given by the following average taken over 
one cyclotron period: 



A(t) 



,ix(*) 



, *(*) = — f dt'8B[v{t')] . (24) 



The integral which determines the random phase x(^) 
should be done on the exact trajectory r(t). Note that 
Eq. (|23| ) gives the nonadiabatic shift both along and 
across the drift trajectory. Since only the latter is of 
interest, one should project the result of the integration 
( p3| ) onto the axis perpendicular to the direction of the 
drift of the outgoing particle. 

Since the nonadiabatic mixing is determined by the 
short wavelength Fourier harmonics of the perturbation 
[Eq. (|23|)1, it is the analytical properties of the function 
A(t) and, therefore, of the correlator (SB(0)5B(r)) that 
are important. In the CF problem, this correlator has 
branch points as a function of r at r = ±2id. However, 
in order to calculate the scattering probability, which is 
given by Eq. (|23|), one has to find the singularities in 
5B[r(t)] as a function of time t and average the result. 
For a given perturbation SB(r) this purely mechanical 
problem of finding the Fourier asymptotics of the inte- 
gral along the path r(t) may be quite complex, but we 
can circumvent the difficulties by performing the con- 
figurational averaging first. As was already mentioned, 
the effective scatterers are characterized by anomalously 
large fluctuations of the drift velocity 



v d (s)=v F (R c /2B)\VB(s)\ , 



(25) 



where s is the coordinate along the path. To see this, 
one can use the exact solution of the single-scattering 
problem considered above. Let us first assume, for the 
purpose of illustration, that the large Vd(s) does not 
change appreciably on the scale of d. Equation then 
tells us that a single impurity located on the trajectory 
which passes through the fluctuation with large V4 yields 



Ap(vd) oc exp(— dw c /vd)- One sees that Ap(vd) grows ex- 
ponentially with increasing Vd- Now, the linear density of 
the fluctuations with large Vd along the percolating path 
is of order p(vd), where the Gaussian probability that the 
drift velocity at a given point is larger than Vd reads 



p(v d ) = cxp{-vj/2(v 2 dx )); 

1 2 \ 3 2 / Rc\ 2 ( b q ^ ' 

Averaging [Ap(vd)] 2 with p{vd), we thus get 



(\Ap{v d )] 2 ) oc exp(-S inin ) ; 



S n 



d 2 L0 2 



2<*4> 



(26) 



1/3 



The "optimum" drift velocity which determines this av- 
erage is v° d = (v 2 dx ) 1/2 (4d 2 w c 7 (v 2 dx )) 1/6 . As is clear, 

S> (Vd x ) at B ^> Bo- The optimum fluctuations 
yield the Gaussian behavior of the scattering rate: 



S mi „ - c{B/B ) 2 , B/B Q » 1 



(27) 



with the coefficient c = I8 1 / 3 ~ 2.62 (here we assumed 
a = 1/ v2) • This simple derivation of the exponential 
dependence of ({Ap(vd)] 2 ) captures the essential physics 
and yields a correct parametric estimate for S m i n ; how- 
ever, it is not exact in that v d {s) in the optimum fluc- 
tuation is not, in fact, constant on the scale of d, and 
for this reason it does not give the correct value of the 
numerical coefficient c in ( p7| ) . To obtain the asymptoti- 
cally exact numerical coefficient in S^in, we have to use 
the optimum-fluctuation method in the whole configura- 
tional space. The optimum configuration is characterized 
by the function v d {s) and the shape of the drift trajec- 
tory. We write the phase factor e lUct in Eq. (|2^) as e I¥,< - s ', 
where we introduce the s dependent phase 



ip(s) = UJ C 



ds' 
Vd(s') 



(28) 



and notice that the exponent 5 m i n is determined by the 
phase ip(s) picked up at the singular point of the pertur- 
bation SB[r(t)] regarded as a function of the longitudinal 
coordinate s. As can be verified by varying the shape of 
the trajectory, the minimum "action" S m [ n is acquired 
along a straight path and the quantity to be calculated 
is therefore 




(29) 



where the integral should be done along the straight line 
connecting the points s = —id and s = id in the com- 
plex plane of the variable s [[30|. This average deter- 
mines, with exponential accuracy, the diffusion coefficient 
D± oc cxp(— 5 m in), which is defined as 
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D 



lim t 



Ap(t)Y 



(30) 



where Ap(t) is the nonadiabatic shift across the perco- 
lating drift trajectory in time t. The exponent can be 
written as a sum of two terms, S m i n = W\ + Wi, where 



and 



W 2 = iu>. 



{VdxVdx) 



2 J (2tt) 2 dx i q dx -i 
ld dx 



e 2 B 1 



(31) 



(32) 



is the Fourier transform of the drift-velocity correlation 
function T(r) = (vd x (0)vdx(r)) ■ Here W\ determines the 
probability for the optimum fluctuation v dx (r) to occur, 
while W2 describes the nonadiabatic scattering on this 
fluctuation. The variational equation SW/Sv dx = yields 



v° dx (r) 



id 



dx' 



T T(x - x',y) , (33) 



which is a non-linear integral equation for v dx (x,0) with 
x G (—id, id). Its solution defines, by means of analytical 
continuation, the optimum fluctuation on the real axis of 
x. Dimensional analysis of this equation shows that the 
solution has the form v dx (r) — vp(B n /B)Q(r/d), where 
Q is a dimcnsionlcss function of order unity, which leads 
again to Eq. (p7|). To find the exact value of c, one has 
to determine the function Q , which requires solving the 



integral equation (33). We used a variational approach, 
choosing the trial function v dx (r) = «r(r) with the vari- 
ational parameter k. This trial function is the optimal 
fluctuation for the (slightly different) problem of finding 
a large value v dx (0) at the point r = and should give 
a good estimate for c. The result is c ~ 2.28, close to 
the value found above within the simplified considera- 
tion neglecting the spatial variation of Vdx on the scale 
of d. 



C. Magnetotransport at a ~ 1: Conductivity of the 
composite fermions at v 7^ 1/2 

We are now prepared to calculate the conductivity at 
B ^> Bq. Similarly to the case of the snake-state perco- 
lation, the nonadiabatic transitions create a conducting 
network with the elementary cell of size £ ~ d(d/Rd) 4 ^ 3 
and perimeter L ~ ^d/ Rd, where Rd, the width of the 
links of the network composed of the drift trajectories, 
obeys the equation 



This equation is the condition of connectivity of the net- 
work. The conductivity due to the nonadiabatic mixing 
of the drift trajectories is thus given by 



? 9 

— ^-VdRd oc —kpd exp 
fir h 
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(35) 



The dissipative transport is seen to be strongly sup- 
pressed beyond the scale B ~ Bq. Let us emphasize 
that the crossover to the adiabatic regime occurs in the 
CF system at a small deviation from half- filling v — 1/2: 
B ~ Bq corresponds to the shift \v — 1/2 1 ~ \jkpd <C 1. 




R 2 d ~ D ± L/v d 



(34) 



co c d/v F 

FIG. 8. Magnetoresistivity at a = 0.2, 0.35, and 0.5; the 
cyclotron frequency in the scaling of the a;-axis refers to the 
external magnetic field, cu c — eB/mc. The full, dashed, and 
dotted lines show the experimental magnetoresistivity around 
v = 1/2 according to the data of Refs. @, and |§, 

respectively (in this case uj c refers to the effective magnetic 
field B = B — 2hcn/e). The sample parameters (carrier den- 
sity n, undoped spacer width d and zero-field mobility p) are: 
n = 1.53 x 10 11 cm" 2 , d = 80 nm, p = 5 x 10 6 cm 2 /Vs 0, 
n = 0.86 x lO^cm" 2 , d = 80 nm, p = 10.5 x 10 6 cm 2 /V s [ 32 1 , 
n = 1.31 x lO^cm -2 , d = 120 nm, p = 3.5 x 10 6 cm 2 /!/ s [26|. 



Calculation of the conductivity tensor components at 
B < Bq cannot be done analytically at a ~ 1. We have 
performed extensive numerical simulations to study the 
overall shape of the magnetoresistivity. The results for 
a = 0.2, 0.35, and 0.5 are shown in Fig. ^ in comparison 
with experimental data on the magnetoresistivity in the 
vicinity of v = 1/2 from Refs. We recall, that 

although a simple model of uncorrelated impurities with 
the concentration rii = n gives a = 1/ the actual 
value of a may be somewhat smaller because of impurity 
correlations. As is seen, the experimental curves are rea- 
sonably close to the numerical data with a ~ 0.25-^0.35. 
The numerical results show a positive magnetoresistance 
at B < Bq, which is followed at larger B by a fall-off 
in agreement with the analytical prediction (|2^). Both 
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these features are clearly observed in the experiment. 
Better agreement between the theory and the experiment 
is hardly possible, taking into account the difference be- 
tween the experimental data obtained by different groups 
and on different (though nominally very close) samples. 
Sufficiently far from half-filling the experimental curves 
start to show magnetooscillations, which have not been 
included in the classical model above. We will discuss 
the issue of the magnetooscillations below. 



D. Magnetotransport at weak disorder, n<l. 

We now turn to consider the case of weak disorder 
a <C 1. The Gaussian behavior of the conductivity at 
large B is in fact a general property of the nonadiabatic 
transport in the RMF. We can use the same optimum 
fluctuation method as for the CF system above to get 



S min = c(2a 2 ) 2 ^ {B/Bo 



B/Bq > a~ 



This result is valid if B is large enough; namely, it re- 
quires that the shift of the guiding center in the field of 
the optimum fluctuation after one cyclotron revolution, 
<5, and the Larmor radius R c be both smaller than d. 
The condition R c <C d breaks down with decreasing B 
at B/Bq r~j a" 1 . Provided a < 1, the exponent j3^ ) 
taken at the crossover point is large, which means that 
5 <C d and we are still deep in the adiabatic regime. Next 
we consider the regime B/Bq < a^ 1 . What changes at 
smaller B is that the drift velocity of the guiding center at 
point r is now determined by an effective RMF 5B cS (r), 
the amplitude of which is smaller than Bq because of 
the averaging of the fluctuations over the large cyclotron 
radius. We define SB (r) by writing the general expres- 
sion for the drift velocity to first order in SB/B in the 
form 

v d (r) = -— J [JB(r + R v ) x R v ] (37) 

= v F ^ \VSB cS {r) x B] , 
2B 

where R v — R c (cos<p,sm(p) and ip is the phase of the 
cyclotron rotation around the point r. Stokes' theorem 
then yields 



5B cS (r) = -j-j / d 2 R8B(r + R) 

n R c Jr<r c 



(38) 



Equation (^) tells us that SB eS is given by the magnetic 
flux through the cyclotron orbit, so that the drift occurs 
along the lines of constant flux (not the lines of constant 
magnetic field averaged along the orbit). If SB(r) is a 
smoothly varying function on the scale of R c , 5B cS (r) 
coincides with SB(r), otherwise the averaging leads to a 
strong suppression of the fluctuations of SB cS (r) . More 



specifically, at R c ^> d the field SB cS (r) is characterized 
by two spatial scales: it has a short-range component 
whose correlation radius remains of order d - its charac- 
teristic amplitude is 



B\ 



eff 



B {d/R c f' 2 



(39) 



- and a long-range component which has a larger ampli- 
tude ~ Bod/R c but fluctuates on the much longer scale 
of R c . Since the drift velocity is given by the gradient 
of 5B cS (r), it is the short-range fluctuations that deter- 
mine V(j(r). Applying the optimum fluctuation method, 
we now have 



a 



1/3 



B/B a 



1/3 < b / Bq < a -i ( 40 ) 



(36) E 



At B/Bq ~ a -1 / 3 the factor e~ Smin becomes of order 
unity and at smaller B the transport is described, in the 
first approximation, in terms of the conventional Drude 
theory; namely, the conductivity grows with decreasing 
B as (7 X x ~ kpd(Bo/B) 2 until it saturates at B/Bq ~ a 
at the level given by Eq. (||). 

Let us briefly discuss the (weak) deviations from the 
Drude behavior at low magnetic fields B. The Drude 
theory predicts zero magnetoresistance, Ap xx /Pxx = 0. 
The low-field magnetoresistance ataCl was studied in 
|uJEl|, where it was found that 



-3a I lnal 



,B 

B 



B < aB< 



o ■ 



(41) 



This result is valid in the weak magnetic field regime 
u c T tr <^ 1 (which is equivalent to the condition B <C 
aBo). As our numerical results indicate, this (very weak) 
negative magnetoresistance crosses over to a much more 
pronounced positive magnetoresistance in the intermedi- 
ate range aB <C B <C a~ 1 / 3 i?o, so that the resistivity 
shows a maximum at B ~ a -1 / 3 £?o before it starts to 
drop exponentially according to Eqs. (35), (36), (|4C 



This positive magnetoresistance remains for all a be- 
low ~ 0.5, see Figs. [|, ^, and is strikingly similar to 
the experimentally observed positive magnetoresistance 
of composite fermions near v = 1/2. We are currently 
working on a theoretical explanation of this effect. 



E. Magnetotransport at strong disorder, a 3> 1: 
Localization of the snake states at B > B c 



Let us address now the finite- B transport in the op- 
posite limit of strong disorder a 3> 1. The asymptotics 
( p6|) remains unmodified at sufficiently large J5, namely, 
at B » Bo/a 1 / 2 , where R c <C R s - At smaller B 7 the tra- 
jectories close to the zero-SB lines become localized in 
the adiabatic traps [Sec. [ij] and start to transform into 
the snake states. A new feature is thus the appearance of 
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a competing mechanism of the conduction - the snake- 
state percolation. Indeed, the snake-state trajectories are 
tied to the lines of B(r) = 0, where B(r) =B + 5B{r) 
is the total field. On the other hand, the percolating tra- 
jectories are those that follow the lines of SB(r) = 0, 
independent of B. At large B ^> Bo/a 1 / 2 , the percolat- 
ing and snake-state trajectories are therefore separated 
in space: the snake states are closed orbits localized deep 
inside the elementary cells of the conducting network. 
Now, at B < So/a 1 / 2 , the nonadiabatic scattering is due 
to two mechanisms: the exponentially weak nonadiabatic 
corrections to the dynamics of the snake-state angle 9 and 
the breakdown of the adiabaticity of the snake states at 
the saddle points. The latter mechanism leads to the for- 
mation of percolation clusters. Thus, there exists a well- 
defined value of B = B c below which the snake states 
form a continuous network and can percolate through 
the entire system. The critical field B c is of the order of 
B ps/d, i.e., 



B r , 



Bna 



-2/3 



(42) 



which is the characteristic amplitude of the magnetic field 
at the critical saddle points. It is worth stressing that 
there is a clear separation of the adiabatic and snake- 
state regimes: the nonadiabatic scattering within the fi- 
nite clusters only leads to an exponentially narrow uncer- 
tainty in the position of the critical point. If one neglects 
this nonadiabatic smearing of the transition, the conduc- 
tivity can be expressed as 



e kpd 



(43) 



where T{x) is a dimensionless function, such that ^"(0) ~ 
1 and T(x > 1) = 0, and o~ xx at zero B coincides with 
that in Eq. (16). The magnetoresistivity for a = 1.5 and 
a — 4 is shown on a logarithmic scale in Fig. ^ (for com- 
pleteness, we also included the data for smaller values of 
a, which have been already displayed in Fig. |J). We see 
that the resistivity indeed shows an exponential fall-off 
beyond a characteristic field consistent with Eq. ([42|). 



To find the critical behavior of o~ xx oc (B c — B) 1 near 
B c , we formulate an auxiliary percolation problem in 
more conventional terms. Consider equipotential con- 
tours in a random potential with a characteristic ampli- 
tude Vo and the correlation length d and pick up all con- 
tours within the energy band (—A, A), where A <C Vo. 
These contours form a percolation network, the size of 
the elementary cell of which is £(A) ~ d(Vo/A) 4 / 3 . The 
characteristic width of the links of the network is of order 
dA/Vo- This is a standard percolation problem. Now, 
let us shift the energy band corresponding to the perco- 
lation network: namely, we introduce a parameter e and 
consider all contours within the band (—A + e, A + e). 
Clearly, the system undergoes a percolation transition at 
e = ±A. In this new percolation problem, the length 



£(A) plays a role of the elementary scale, so that the 
characteristic radius of the percolation cluster is now 
£(A,e) ~ £(A,0)(A/|A - e|) 4 / 3 (for positive e). Next, 
let particles propagate along the equipotential contours 
"ballistically" . As we know already from Sees. 0g, 
the conductivity a of such a network scales as the typ- 
ical width of the conducting links d(A — e)/Vo, i.e., 
o oc (A — e). Assuming that our original problem can be 
mapped onto the problem above, we will get the critical 
exponent for the conductivity of the snake-state network 
t = 1, i.e., T(x) ~ 1 — x at x — ► 1. 
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FIG. 9. Magnetoresistivity at a = 0.2, 0.35, 0.5, 1.5, and 
4. The dotted lines show the critical field _B e , at which the 
exponential fall-off begins for a > 1; see Eq. (143). 



F. Magnetooscillations 

Until now we have dealt with the transport properties 
of classical particles. The classical description is justi- 
fied by the large value of the parameter kpd 3> 1: first, 
this parameter enabled us to calculate the conductivity 
by examining the microscopic dynamics in terms of clas- 
sical trajectories; second, it guaranteed the existence of a 
wide range of the fields Bq and B in which the thus cal- 
culated conductivity a xx ^ e 2 /h [cf. Eqs. It is 
the latter condition that allowed us to neglect the quan- 
tum interference of multiply scattered waves and related 
localization effects. In this section, we discuss the small 
quantum oscillations of function of B, which 
are conventionally termed the Shubnikov-de Haas (SdH) 
effect. We argue that the physical picture of the magne- 
tooscillations in the smooth random magnetic field (and, 
in fact, in the limit of a long ranged random potential 
as well) is rather peculiar from the conventional point of 
view. 

We first recall the standard results for the case of a ran- 
dom scalar potential with a sufficiently short correlation 
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length. Within the usual approach |33|.[54]], the magne- 
tooscillations of a xx are related to those of the total den- 
sity of states pf- o-° s x /a xx oc pp c /pF, where a° s x and p^ sc 
denote the oscillating parts of o~ xx and pf- The damping 
of the oscillations in low fields is described in terms of the 
single-particle relaxation time r s (also termed sometimes 
the quantum relaxation time) ]3^- |3^ , p^ | , which is equal 
to the transport mean free time r tr for the white-noise 
random potential and is ~ r tr / '(kpd) 2 for the random 
potential with correlation length d > kp l . Specifically, 
the exponent of the Dingle factor e~ s is S = tt/uj c t s , 
so that at lo c t s ~ 1 the oscillations of both the density 
of states at the Fermi level and the conductivity become 
strong. We note that this crossover to the strong os- 
cillations occurs in the case of short-range disorder at 
&xx/(e 2 /h) ~ kFl(T s /Tt r ) 2 3> 1, assuming that the ratio 
r fr /r s is not too large. With further increasing B, at 
uj c t s ^> 1, the density of states exhibits a series of peaks, 
well separated from each other, so that in the centers of 
the valleys between the peaks the quantum localization 
starts to develop rapidly, which leads in turn to the ap- 
pearance of quantum plateaus in the Hall conductivity. 
The full Hall quantization takes place, however, at much 
larger B, when the value of the conductivity in the cen- 
ter of the peak drops to a value ~ e 2 /h. Therefore, in 
the short-range random potential there exists a paramet- 
rically broad region between the appearance of the SdH 
oscillations and the fully developed quantum Hall effect. 

In the case of a white-noise random magnetic field, 
the situation is different. The short-range fluctuations of 
the magnetic field are accompanied by long-range fluc- 
tuations of the random vector-potential. As a result, 
the single-particle relaxation rate now diverges due to 
the strong small-angle scattering |^,|l9|. The divergence 
is cut off by the characteristic length scale in the prob- 
lem, which is the cyclotron radius. As a consequence, the 
damping factor S takes a different form, S = AttEf /u 2 Tt r 
fL9f . It follows that at S ~ 1, where the oscillations 
become observable, a xx ~ e 2 /h. This is, therefore, a 
marginal case: the SdH oscillations become appreciable 
in the same region of B, where the quantum localization 
effects get strong. 

Now let us consider the limit of the long-range RMF. 
In this case, the oscillations of the total density of states 
are damped exponentially strongly in the whole region, 
where the classical conductivity is > e 2 //i, and are thus 
of no importance. However, the crucial thing to notice 
is that the oscillations of a xx and of the total density 
of states pf are no longer directly related to each other. 
Indeed, as shown above, a xx is determined with growing 
B by a progressively smaller fraction of the total number 
of trajectories. In the adiabatic limit, only trajectories 
close to the zero-<5-B contours contribute to a xx . It follows 
that the oscillations of associated with the oscil- 

lations of the density of these conducting states only. In 
particular, this means that the SdH effect and the magne- 



tooscillations of thermodynamic quantities (de Haas- van 
Alfven effect) will be characterized by completely differ- 
ent damping factors. 

We now turn to the Dingle factor for o~ xx . Quite gen- 
erally, the SdH effect is due to the quantum interfer- 
ence of two waves propagating along quasiclassical tra- 
jectories for which the number of cyclotron revolutions 
is different by ±N for the Nth harmonic of the oscil- 
lations Thc Dm s le factor for the first ( and 

most prominent) harmonic with N = 1 can therefore be 
represented as 

e~ s cos^ = G a cos$ Q ^ , (44) 

where a labels trajectories in a given realization of dis- 
order, G a is the weight with which the trajectory a con- 
tributes to a xx , J2 a G a = 1, <§> a is the phase which is 
acquired by a particle moving along the trajectory after 
one cyclotron revolution, -0 is the phase of the SdH oscil- 
lations, and ( ) denotes ensemble averaging. In general, 
the phase factor cos $ should be averaged both over dif- 
ferent trajectories and along one trajectory. However, 
in the limit of smooth disorder, the action $ a , which is 
the dimensionless magnetic flux through the cyclotron 
orbit, is the adiabatic invariant characterizing the tra- 
jectory. A subtle point here is that the van Alfven drift 
occurs along the lines of constant flux <!> - not the lines 
of constant field B [see Eq. (psj)]. At first glance, this 
difference might seem to be irrelevant in the case of long- 
range RMF with d^> R c . In fact, however, it is of crucial 
importance for the calculation of the amplitude of the os- 
cillations. Let us illustrate this point by first deriving S 
along the following line of argument. We know already 
that the extended trajectories form a percolating net- 
work along the contours of zero SB(r). Since the width 
of the links of the conducting network is much smaller 
(in fact, as shown above, exponentially smaller) than R c , 
it is a good approximation to place the guiding center 
on the contour of zero SB(r), calculate the RMF flux 
through the cyclotron orbit, and average e I<5 * over 
different positions of the guiding center on the zero-(5-B 
line. This would give a contribution to the Dingle-factor 
exponent S = ((6<I>) 2 ) /2 ~ (eB Q /hc) 2 {R 2 /d) i (Tj]. In 
fact, however, the actual trajectory of the guiding center 
is slightly shifted from the contour of zero 6B(r) - by 
an amount which exactly cancels the above contribution 
to S - since it is the flux that is the adiabatic invariant. 
It follows that the average over the flux should be done 
only over different trajectories: $ a in Eq. (|44|) is, in the 
adiabatic limit, constant for a given a. 

Having established the conservation of the flux $ a 
along the trajectory we can rewrite Eq. (|44|) in the fol- 
lowing form 

e~ s cosip = / d$G(>I>) cos $ , (45) 
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where G(<1>) = (J2 a G a 6(® - $ Q )) is understood as the 
ensemble average taken across a link of the conducting 
network. The function g(&) is represented as a narrow 
peak of width A$ centered at $o = 2ir(mvp /2Tiu> c ). 
The dimensionless magnetic flux A$ can be expressed 
in terms of the width Rd [Eq. (pif)] of the link of the con- 
ducting network by A$ ~ (e/hc) AB R 2 , where AB ~ 
B^Rd/d is the characteristic change of the RMF across 
the link. We thus see that the broadening A$ is related 
to the conductivity of the network [Eq. ([35])] and obeys 
the equation 



cr_, - — A$ 



(46) 



This rather remarkable result implies that the Dingle 
factor e~ s is a function of the single variable g = 
a xx /{e 2 /h): 



dxQ(x) cos(27rgx) . 

-oo 



(47) 



Here, we have expressed G($) = (27r<7) _1 <3[(<I>— $o)/27r0] 
in terms of the parameterless function Q(x) which falls off 
at \x\ ~ 1 and is normalized according to J dxQ(x) = 1. 
According to Eq. (^7|), the Dingle factor is represented 
as the Fourier transform of the smooth function Q(x). 



It is worth noting that Eq. (47) can be interpreted also 



in terms of the local Landau levels. In the language of 
the quasiclassical quantization, the contribution Gpj of 
the Nth Landau level to the conductivity a xx = ^ N Gn 
falls off beyond the band of width AiV ~ g ^> 1 around 
N = mvp/2huj c , where the number of effectively con- 
ducting Landau levels is determined by the change of the 
flux across the link ; AN = A$/2n. With this fact taken 
into account, Eq. ( f46| ) takes the familiar form: the con- 
ductivity is of the order of e 2 /h times the number of the 
conducting channels in the effective network. Applying 
Poisson's formula to the sum Gpj we again arrive at 
Eq.©. 

Hence, the SdH oscillations due to the oscillations of 
the density of states of the "conducting" particles be- 
come observable at g ~ 1. However, as has been already 
mentioned, there exists another effect, which leads to the 
appearance of the magnetooscillations, — the quantum 
localization. According to the scaling theory of the quan- 
tum Hall effect 0, 



Sqhe(s) = 2ng 



g > 1 



(48) 



irrespective of any microscopic details, in particular, the 
value of a. Equations J47|), (Esj) tell us that both types of 
oscillations become observable at g ~ 1. To decide which 
oscillations are stronger, one should calculate the Fourier 
asymptotics (47). This requires knowing the precise 
shape of the function Q{x). It appears that in the case 
of the percolation network Q(x) can be obtained only by 



means of a numerical simulation. Here, we restrict our- 
selves to concluding that the number of oscillations ob- 
served scales in the CF problem as p c ~ kpdj \v}^ 2 (kFd). 
Since g oc fc^ciexp[— (p c /p) 2 ], where p = Ep/huj c , the os- 
cillations disappear extremely fast with increasing p (i.e. 
decreasing B). These findings are in agreement with ex- 
perimental observations. The typical number of oscil- 
lations observed in the best samples with kpd ~ 15 is 
p c = 7 -r 9. The oscillations start indeed to develop at 
g i=a 1, as shown in Fig. [l^, where the experimental data 
of Ref . |^6| are represented in terms of the CF conductiv- 
ity. Finally, the damping of oscillations with decreasing 
B is extremely fast, so that the Dingle plot is strongly 
non-linear pa,E6[. 
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FIG. 10. Composite fermion conductivity a xx as a func- 
tion of the effective magnetic field extracted from the data of 
Ref. ^|. The magnetooscillations start to develop when the 
conductivity (in units of e 2 /h) drops down to a value ~ 1. 



IV. AC TRANSPORT 

The frequency-dependent dissipative conductivity in- 
volves the Fourier transform of the retarded velocity- 
velocity correlator: <j xx (uj) — e 2 pF (\v x \ 2 J ), where 



(k| 2 ) =Re / dte 1 ^ (v x (Q)v x (t)) 



(49) 



In the quasiclassical limit, which we consider here, the 
integral over t is understood as the integral along a clas- 
sical trajectory, characterized by the velocity v(t), while 
( ) denotes the averaging over the trajectories. The aver- 
age is taken over all trajectories of electrons at the Fermi 
level. However, as discussed above, if the magnetic field, 
either B or Bo, is strong enough, most trajectories do 
not make a significant contribution to the dc conductiv- 
ity and the transport is governed by the percolation of a 
small portion of their total number. In this section, we 
turn to the ac conductivity. We will show that changing 
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the frequency, one effectively probes the motion on the 
percolating network on different spatial scales. This re- 
sults in a strong non-Drude frequency dispersion of the 
conductivity. 



A. ac conductivity at a ~ 1, B > Bq: application to 
the composite fermions around v = 1/2. 

We start with the ac conductivity at a ~ 1 in a strong 
external field B ^> Bq (which is relevant to the CF 
problem at sufficiently strong deviation from v = 1/2, 
Av 1/kFd). Our analysis in Sec. [II relies on the 
characterization of the electron dynamics at B 3> Bq by 
three degrees of freedom with different frequency scales: 
the fast cyclotron rotation, the slow drift, and the still 
slower nonadiabatic diffusion across the drift trajecto- 
ries. It is the slowest degree of freedom, the nonadiabatic 
diffusion, that yields unbounded (extended) paths. Let 
us first consider the contribution to the velocity-velocity 
correlator (49) which comes from these extended trajec- 
tories, (v x (0)v x (t)) cxt . Following the approach developed 



in Sec. we parametrize it as a function of two vari- 
ables 



(v x (Q)v x (t)) ext =v 2 d -F v ( 



v d t v d t 



), 



(50) 



where v d = (3/8) 1 ' 2 VFR c Bo/dB is the r.m.s. drift ve- 
locity, £ oc and L oc R^" 1 are the characteristic 
size and perimeter of the elementary cell of the perco- 
lation network, and v = 4/3. The factor d/£ gives the 
partial density of the percolating states (i.e., the por- 
tion of the area £ x £ occupied by the trajectories of 
size f, LR d ji 2 ~ d/f). At d/v d < t < L/v d , the scal- 
ing function exhibits a power-law behavior F v (t,0) ~ 
t~ x reflecting the fractal dimensionality of the links of 
the network. The exponent x can be found by equat- 
ing (v x {Q)v x (t)) 1 1 and the effective diffusion coefficient 
£ 2 (t)/t, where £(*) ~ d[L(t)/d] u ^ v+ ^ and L(t) ~ v d t, 
which yields x = 2/{v + 1) = 6/7. The network model 
we have used is justified by the condition x < 1 (the 
Harris criterion, v > 1), so that the integral over t in Eq. 
(kM) at ui = is determined by t ~ L/v d . We thus have 



(v x (0)v x {t)) ext ~ v, 



,d 

'£ \vdt 



-<t<- 

v d ~ ~ v d 



(51) 



Now, to calculate the frequency-dependent correction to 
the dc conductivity at u> < Vd/L, we need to know the 
behavior of the correlator ([50]) at t ^ L/v d . If the dif- 
fusion over the percolation network were completely un- 
corrected, F v would decay exponentially with increasing 
t. We will argue below that in fact F v changes sign at 
t ~ L/vd and falls off as a power law at larger t: 



F v (r,r')^-(r'/rr/T' 2 



T > t' > 1. 



For the correlator d5G) this gives 



(v x (0)v x (t)) ( 



d£ 
t 2 



t> 



v d 



(52) 



This long-time tail in the correlation function is similar 
to the one found in the Lorentz gas (noninteracting clas- 
sical particles scattered by a random array of hard discs) 
p7[ . We will study the long-time correlations microscop- 
ically for the realistic case of a weak long-range potential 
elsewhere; here, we introduce a simple phenomenologi- 
cal model suitable for the qualitative description of the 
percolation network. 

Consider the diffusion equation with an inhomoge- 
neous diffusion coefficient D(r). The diffusion current 
is then given by Z3(r)Vn(r, i), where n is the concentra- 
tion of particles. The diffusion propagator n wq in (lj, q) 
space obeys the equation 



d 2 q' 



(qq')Ai-q'' 



1 



(53) 



We now write D(r) = D + SD(r), where D is the mean 
value of D(r), expand the propagator in powers of SD, 
and average Eq. (|3|) over the fluctuations with the cor- 
relator (5D(0)SD(r)) = KD 2 g(r), where g{r) is a dimen- 
sionless function of order unity which falls off on the scale 
of do. The use of the diffusion equation implies that the 
correlation radius do 3> I, where I is the mean free path. 
To first order in A we have for the correction to the av- 
erage diffusion propagator 



d 2 q' 



Sq-q' (qqTZW , (54) 



where 2? wq = (— ilu + Dq 2 )^ 1 . From Eq. ( |54| ) we de- 
duce the uj dependent correction to the conductivity 
Sa xx (uj) = e 2 ptu 2 lim^o q~ 2 ReSn U q: 



Sa xx {uj) 



d 2 q 

(2nf 



3q 



Dq 2 



-iu> + Dq 2 



(55) 



i.e., 5a xx {us) / o xx ~ kg$\uj\/lQD at small w [we drop 
here an uj independent term in 5a xx (u>)], which implies 
a t~ 2 long-time tail in the velocity-velocity correlator in 
Eq. H). 

On the percolation network, the size of the effective 
scatterers and the effective mean free path are both of 
the order of £, i.e. do ~ I ~ £. Furthermore, strong fluc- 
tuations of the geometry of the percolating cluster imply 
that A ~ 1. We expect that the 1/t 2 tail, the existence of 
which has been demonstrated above in the phenomeno- 
logical model with do ^> I and A<1, will not disappear 
if we set d ~ I and A ~ 1. Substituting these estimates 
into Eq. (|55|) we get 



8a xx (uj) \uj\L 



< 



I'd 



Vd 

L ' 



(56) 
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which corresponds to Eq. (52) ps[| . We recall that the 
characteristic scale L entering Eq. (B6ft is given by 



L~ d 



me 2 Vdd 

®~XX 



7/3 



\ </exp ( —S min 



(57) 



It is worth noting that the classical kinetic correla- 
tions compete with the quantum ones and win, unless 
the frequency is exponentially small. Specifically, as is 
well known, the quantum localization effects in 2d lead 
to a t~ x tail in the correlator ( pf9| ) and, correspondingly, 
to a In correction to the conductivity. This quantum 
correction is of special interest because of its divergence 
for u> — ► 0, in the thermodynamic limit. The classical 
correction, proportional to |cj|, does not diverge, but it 
is also interesting, both theoretically and experimentally, 
since it is nonanalytical in u and is much larger than the 
quantum one even at very low to. The point is that the 
localization correction is a series in powers of the small 
parameter 1/kpl, where kp is the Fermi wave vector, 
while the relevant parameter for the classical corrections 
is d/l, where d is the correlation radius of disorder. If the 
disorder is long-ranged (kpd ^> 1), the classical correc- 
tions are dominant in a wide range of to. We will discuss 
the classical corrections in more detail elsewhere. 

Note that a xx (ui) behaves nonanalytically at half- filling 
also in the integer quantum Hall regime. At the integer 
quantum Hall transition, the frequency dependent cor- 
rection 5a xx ~ | y / 2 is related to corrections to scaling. 
The leading irrelevant scaling exponent y at the quan- 
tum phase transition was found numerically to be equal 
to 0.38 ± 0.04 and 0.35 ± 0.05 in Refs. @ and @, re- 
spectively. In |42| y was argued to be equal to 77, where 
77 ~ 0.4 ^,^3|is the critical exponent of eigenfunction 
correlations. However, in a long-range random potential 
there exists a classical non-analytic term Sa xx ~ |w| yc1 / 2 
with y c \ w 2 [fill , which dominates 8u xx {uj) in a wide 
range of frequencies. 

Coming back to Eq. (|56|), we see that the frequency- 
dependent correction becomes strong at oj ~ Vd/L. We 
now turn to higher frequencies. The scaling form (51 ) im- 
plies that the contribution to <j xx {u>) from the extended 
trajectories behaves as lu~ 1+x at lu ^> Vd/L, i.e., it slowly 
decreases as w -1 / 7 with increasing frequency. Let us 
show that the extended trajectories do not determine 
the conductivity any more and the main contribution to 
&xx{u) comes now from localized drift trajectories. We 
will show that in fact <j xx {uj) grows with to. We neglect 
the weak nonadiabatic scattering between the closed drift 
trajectories and represent the velocity-velocity correlator 
(p9|) for the periodic drift orbits in the form 



(Ml) 



loc 



f T dt , , 

/„ T Z{t) ' 



27m s 



Here r(t) is a closed trajectory with the period T and 
the angular brackets denote the average over both T and 
the shape of the trajectory at given T. The average is 
determined by trajectories with T ~ u>~ , their perime- 
ter and size are L w ~ Vd/to and £ w ~ d{L u j 'd) v '( u+1 >, 
respectively. The partial density of states corresponding 
to these trajectories is ~ pd/^ w [see the paragraph after 
Eq. (pQ)]. The estimate for the correlator (p8) thus reads 



<k£) loc ~w 2 £xd/6,x M- 1 



(59) 



where the factor d/£ u stands for the partial density of 
states and the last factor w^ 1 comes from the averaging 
of the delta functions in Eq. (|5"s|). This yields the ac 
conductivity of the form 



~ e 2 pv d R u 



(60) 



where R u ~ G?(d/£ w ) 1 / l/ oc |w| 3 / 7 is the characteristic 
width of the links of the conducting network. We get 



®xx 



3/7 



< 



< v -i. (6i) 



Equations (|5^) , ( |6l| ) describe the behavior of the con- 
ductivity at u> ^ Vd/ d. At still higher frequencies, a xx {ui) 
is determined by the velocity-velocity correlations in the 
crossover region between the "ballistic" drift on the spa- 
tial scale much smaller than d and the "diffusive" mo- 
tion over the fractal network on larger scales. At the 
crossover, the velocity-velocity correlator can be param- 
eterized as (v x (0)v x (t)) = v 2 F v (vdt/d), where F v (t) is a 
dimensionless function of order unity [cf. Eq. (^0|)] . Note 
that it is not sufficient at w > Vd/d to know the behav- 
ior of the correlator at t ~ lo^ 1 . Since the conductivity 
is expressed in terms of the high-frequency Fourier com- 
ponent of (v x (0)v x (t)}, which is an even function of t 1 
a xx will fall off exponentially with increasing uj (until 
u> reaches the low-frequency wing of the cyclotron res- 
onance). Therefore, to get the asymptotic behavior of 
<r ra (w), we need to know the analytical properties of the 
correlator as a function of t, i.e., the exact shape of the 
function F v (t), which requires a numerical simulation. 
We expect, however, that the function F v (r) has a simple 
analytical structure with singular points at Imr ~ ±1, 
which yields the exponential falloff of the form 



In a q 



i\d/v d . 



(62) 



(58) 



This exponential decay of a xx is limited from the side 
of large frequencies by the disorder-broadened cyclotron 
resonance which dominates <j xx (lo) at ui ~ lu c = eB/mc. 

The arguments of the last paragraph concerning the 
exponential falloff at large ui are also applicable to the ac 
conductivity in zero (or low) B, with the only substitu- 
tion of the Fermi velocity vp for the drift velocity Vd- In 
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contrast to the Drude (white-noise disorder) case, where 
the velocity-velocity correlation function has a cusp at 
t = leading to the slow 1/lu 2 decrease of a xx (uj), in the 
case of smooth disorder (v x (0)v x (t)) is an analytic func- 
tion of t at t = 0, which implies an exponential decay of 
u xx (uj) at uj ^> vp/d. 




0.0 0.5 1.0 1.5 2.0 

cod/v F 



FIG. 11. ac conductivity at a = 0.35 and for different val- 
ues of the magnetic field B. The number near each curve 
indicates the value of the cyclotron frequency ui c — eB/mc 
in units of vf/cL; to convert these values into B/Bo, one has 
to multiply them by 1/q. The result of Drude theory is also 
shown (dotted lines). The statistical noise in the data at 
uicd/vp = 1.22 is due to fluctuations of the local cyclotron 
frequency. Inset: the low-B data (u} c d/vF = 0.01) on a loga- 
rithmic scale. 



the ac conductivity is completely different: it increases 
non-analytically at small u> in agreement with Eq. (|6l|), 
see Fig. [jj], and then drops exponentially in a higher fre- 
quency range in agreement with Eq. (]6^), see Fig. [l^; at 
still higher frequencies, the cyclotron resonance (smeared 
by disorder) is observed. It becomes difficult to resolve 
reliably the leading non-analytic correction Sa xx oc \lu\ at 
uj — *■ in high magnetic fields, since it is shifted to the 
very low uj range and is masked by the statistical noise 
present in the numerical data. Note that in Figs. |l2| and 
|i~3| we used the frequency scale v^/d, which is a natu- 
ral scale in the regime of high magnetic field. The value 
of a in Fig. |ll| is approximately the one appropriate for 
the CF system. Therefore, this figure represents our pre- 
diction for the ac conductivity of the CFs at half-filling 
(B = 0) and away from it (B ^ 0). 




FIG. 12. Low-frequency behavior of the ac conductiv- 
ity 8a xx (u) — a xx (co) — a xx (0) in strong magnetic field 
for a = 4.04, B/Bo = 1.00 (which corresponds to 
u! c d/vF — 4.04). The straight line corresponds to the the- 
oretical prediction Sa xx oc |w| 3,/7 , see Eq. (|6l|). The same 
non-analytic behavior of a xx (ui) in the Iow-oj range in strong 
magnetic field B takes place for small a and is in particular 
seen in Fig. [L3] for a — 0.2. 



Figure |ll| shows the results of the numerical calcula- 
tion of the ac conductivity for a — 0.35. Significant 
deviations from the Drude theory fit are seen, which be- 
come stronger with increasing external magnetic field B. 
For zero (or low) B the results are still relatively close 
to the Drude theory, except in the tail (for uj ^S> vf/o 1 ), 
where the conductivity starts to drop exponentially (see 
inset), in qualitative agreement with the theoretical ex- 
pectation [Eq. (p3)]. In the intermediate fields B (see 
the curves corresponding to uj c djvp = 0.2 and 0.37) 
the non-analytic dip ( p6| ) around uj = gets clearly ob- 
served. Finally, in a large magnetic field, the shape of 
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FIG. 13. ac conductivity at a = 0.2 and strong magnetic 
field, B/Bo = 4.38 (which corresponds to uj c d/vF = 0.88), 
shows the exponential fall-off ( |62j ) in the intermediate 
range of frequencies. The dotted line corresponds to 
hicr xx = ~2.2\u\d/vd- 



B. ac conductivity at a > 1, B = 0. 

Let us now consider the frequency dispersion of the 
snake-state percolation (a ^> 1, B = 0). Since the mech- 
anism of the percolation was of no importance in the 
derivation of Eq. , we get 



Sa xx (uj) \u)\L s 



vf 



< 



(63) 



(we substituted here vf for Vd as the effective drift veloc- 
ity of the snake states and L s ~ da 14 / 9 for the perimeter 
of the elementary cell of the conducting network) . The 
correction becomes strong at |w| ~ vp/L s . At larger 
to, the conductivity is determined by the snake states 
that are bounded to the closed zero-5 contours of size 
much smaller than £ s . Specifically, the main contribu- 
tion to a xx (tu) now comes from the trajectories of the 
length L u ~ vf/uj. We estimate their contribution at 
M }t vf / L s as 



pR 8 L^9 2 



So, 



x Z'u ~ a xx (0) 



(64) 



vf < , , < «FP 5 



d 2 ' 



where ^ - d{L w / d) v ^ v+1 \ The first factor gives the 
density of states for the snake-state orbits with the pe- 
riods larger than while the last factor is the ef- 
fective diffusion coefficient on the time scale of order 
u> . We see that u> cancels out in this expression and 
a xx (uj) turns out to be of order a xx (0) [Eq. (filf)]. We 
put \uj\ ~ vpps/d 2 ~ VF/da 2 / 3, as the upper limit for the 



frequency range where Eq. (|64|) is valid: strictly speak- 
ing, the decrease of L u at higher \u>\ is accompanied with 
a growth of the phase volume of the snake states that 
participate in the ac transport by escaping the adiabatic 
traps. This means that the angle 9 which should be sub- 
stituted for Q s , Eq. (H), now increases with \uj\ as Cu , 
where ~ ln 1//4 (u^ It follows that 

ct xx (uj) ~ <r xx (0)— , < |w| < ^ , (65) 

i.e., the conductivity grows with increasing ui (but ex- 
tremely slowly), until |w| becomes of the order of vf /d. 

Note that, in the derivation of Eqs. (|64[) , (|65|) , we re- 
stricted ourselves to the snake states and were not con- 
cerned about the contribution of the drift trajectories, 
which requires comment. Indeed, one might think that 
the arguments that led to the power-law behavior of 
a xx {uj) in Eq. ( |6l"| ) could be used here as well. However, 
there is a new feature that makes the percolation at large 
a and zero B distinct from that at a ~ 1 and large B. 
Namely, now there is no characteristic drift velocity Vd 
which is the same for all trajectories. Specifically, Vd at 
large a depends on the typical distance 1Z between the 
drift trajectory and the zero-B contour. We parametrize 
Vd(H) as 

v d (K) = v F F s (J^j , (66) 

where F s (0) ~ 1 and F s (x) ~ —x~ 2 at x » 1. The slow- 
ing down of the drift with increasing 1Z is related to the 
linear growth of the magnetic field as one moves away 
from the line B = 0. Equation (|6(]) tells us that the 
contribution to the ac conductivity of trajectories sepa- 
rated by the distance 1Z from the zero-SB lines scales as 
Vd(lZ)lZ oc IZ^ 1 , i.e., it decreases with TZ. It follows that 
the drift orbits that surround the snake-state trajectories 
do not contribute to a xx (u>) even at large uj Eq]. 

Thus, the overall picture in the snake-state percolation 
regime is as follows: a xx (uj) exhibits a narrow dip around 
u> = 0, increasing linearly with growing u>; this increase 
is saturated at w ~ vp/L s where the w-dependent cor- 
rection becomes strong; finally, on the (parametrically 
larger) scale of u> ~ vp/d, a xx (uj) starts to fall off ex- 
ponentially The latter regime is similar to that at 
large B [see the paragraph after Eq. (O)]. 

We finally comment on the case a ~ 1, B = relevant 
to the CF problem. In this case L s ~ d, so that the range 
of applicability of Eqs. (UJ), ( |65| ) shrinks away and the 
ac conductivity becomes a function of the single variable 
ujd/vF: 



_ . tod 

<J XX {U>) ~ <Txx\y)F u — 
VF 



(67) 



where F u (x) - F u (0) ~ \x\ at \x\ < 1 and F u (x) falls 
off exponentially at larger x (\nF u ~ — |x|). We have al- 
ready discussed the exponential behavior at large u> in the 
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considered regime a ~ 1, B = in the end of Sec. IV A; 
it is clearly seen in Figs. [ll], [li]. As to the nonanalytic 
dip at small lu, our numerical simulations indicate that 
it is almost unobservable at the values of a describing 
the CF conductivity (a < 0.7) at B = 0. Apparently 
the corresponding numerical coefficient gets very small 
for such values of a. With increase of either a (Fig. |l4| ) 
or B (Fig. [ll]) the non-analytic structure gets resolved 
very clearly. 

In Fig. |lj we show the ac conductivity at relatively 
large a = 4.04 in comparison with that at a = 0.5. The 
a = 4.04 curve is strikingly different from the Drudc 
behavior and shows the features discussed above: the 
non-analytic increase in the small lu region, followed by 
a rapid decay consistent with the theoretically predicted 
exponential falloff. At still higher frequencies a broad dis- 
tribution of cyclotron resonances in the local (random) 
magnetic field is observed. A hump at lucI/vf ~ 1.7 
marks the onset of this regime and corresponds to the 
characteristic snake state frequency to ~ ^^(vp/d). To 
see the exponential falloff in a broader frequency range 
more clearly one would have to consider larger values of 
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FIG. 14. ac conductivity in zero B at a — 0.5 and a = 4.04. 
The scaling of the y axis is different for the two curves, as in- 
dicated by the arrows. Inset: non-analytic dip around lu = 
[Eq. (Eh] at a = 4.04. 



V. CONCLUSIONS 

We have presented a detailed analytical and numerical 
study of the conductivity of a 2D fermion gas in a smooth 
random magnetic field, in the whole range of the parame- 
ters a (strength of the random field) , B (mean magnetic 
field), and lu (frequency). While special emphasis has 
been put on the application of our results to the com- 
posite fermion description of the half-filled Landau level, 
they may be equally relevant to the electron transport in 



a real random magnetic field. Below, we summarize the 
main findings: 

1. At zero magnetic field B, the dc transport has a 
totally different character in the regimes of weak 
(a <C 1) and strong (a ^S> 1) disorder. While in the 
former case a xx cx 1/a 2 [Eq. (En] is correctly given 
by the Born approximation, in the latter the con- 
ductivity is determined by the percolation of snake 
states yielding a xx cx 1/a 1 / 2 (up to a negligibly 
weak logarithmic correction), see Eq. (|l6|). Numer- 
ical simulations confirm these analytical findings 
and allow us to find a xx in the crossover region 
a ~ 1, see Fig. [jj. 

2. In strong mean magnetic field B the particle motion 

takes the form of an adiabatic drift of the cyclotron 

orbits. A non-zero value of <j xx in this regime is 

entirely due to exponentially weak non-adiabatic 

scattering processes. As a consequence, the con- 

— 2 

ductivity falls off exponentially, — In a xx cx B , see 
Eqs. (§7|), (|H), and (||). [At a < 1, an interme- 
diate regime appears, — lner^ cx B, see Eq. (^0|).] 
The numerical simulations have allowed us to find 
the shape of the magnetoresistance in a wide range 
of B for the values of a ranging from 0.2 to 5, 
see Fig. |[ The magnetoresistance p xx (B) shows 
a sharp falloff at large B in agreement with the 
analytical results. Furthermore, at a ^ 0.5 this 
falloff is preceded by a positive magnetoresistance 
in the intermediate range of B. The whole shape of 
Pxx(B) at a ~ 0.35 (as well as its absolute value) 
is surprisingly similar to the experimental magne- 
toresistivity in the fractional quantum Hall effect 
around v =1/2 (Fig. ||). 

3. In contrast to the case of a short-range random 
potential, the quantum magnetooscillations of the 
conductivity start to develop in the range of the 
magnetic field B where the dimensionless conduc- 
tivity g = <y xx /(e 2 /h) drops to a value ~ 1, in 
agreement with experiment (Fig. |l^). These oscil- 
lations are not related to those of the total density 
of states (which are damped much more strongly), 
but are determined by the oscillations of the density 
of states of the particles moving on the percolating 
network, as well as by the quantum localization ef- 
fects. 

4. The ac conductivity also shows distinct features re- 
lated to the deviations of the particle kinetics from 
the behavior following from the Boltzmann equa- 
tion. While at a < 1 and B = the ac conduc- 
tivity ct xx (lu) is relatively close to the Drude form 
(except in the tail, where it drops exponentially), 
at large a and/or B the shape of cr xx (u>) becomes 
totally different. Specifically, it shows nonanalytic 
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behavior in the low-frequency range, see Eqs. (56), 
©, (H) and Figs. related to the long-time 

tails in the velocity-velocity correlation function 
and reflecting the strongly non-Boltzmann char- 
acter of the transport in the percolating regime. 
At higher frequencies, <j xx (uj) starts to drop expo- 
nentially (which reflects the "ballistic" dynamics of 
the snake states or drifting orbits on short spatial 
scales), until it reaches the low- frequency wing of 
the disorder-broadened cyclotron resonance peak, 
see Eq. ([32|) and Fig. [l3|. The ac conductivity at 
a = 0.35 (which is in the range of a relevant to the 
composite fermion problem) and different values of 
the mean magnetic field B shown in Fig. [ll] clearly 
demonstrates the anomalies which we expect to be 
observed in the ac transport around v = 1/2. 
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